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Abstract. We prove that the solutions to fast diffusion stochastic porous 
media equations have finite time extinction with strictly positive probability. 



1. Introduction 
Consider the stochastic porous media equation 

(1) 

dX{t) ~ pA(|X|"(i) sign X{t))dt - A{^{X{t))dt = a{X{t))dW{t), in (0, oo) x O, 
X = on (0, oo) X aO, X{0,x)^x on O, 

where p > 0, a E (0,1),^ is a continuous monotonically non decreasing function of 
hnear growth and a{X)dW — X^fcLi l^k^^kd-Pk, t > 0, where {Pk} is a sequence 
of independent real Brownian motions on a filtered probabihty space ($7, {^t}, P) 
and {cfe} is an orthormal basis in L^{0) which for convenience wiU be taken as the 
eigenfunction system for the Laplace operator with Dirichlet boundary conditions, 
i.e., — Aefc — XkCk in O, Ck — on dO, where O is an open and bounded subset of 
R'^, with smooth boundary dO. We shall assume that J^'kLi f^k^k < Equation 
(H)) for < a < 1 is relevant in the mathematical modelling of the dynamics of 
an ideal gas in a porous medium and, in particular, in a plasma fast diffusion 
model (for a = 1/2) (see e.g. 14 ). The existence and uniqueness of a strong 
solution in the sense to be defined below was studied in [I] , [2] , [3] , [S] for more general 
nonlinear stochastic equations of the form In [3] (see also [I]) it was also 

proven that for a = and d = 1 the solution X = X{t,x) to ^ has the finite 
extinction property: P(r < n) > I — (JJ' e^'-^^^ds) for |a::|_i < C']^^frf 

where t = inf{i > : \X{t,x)\-i = 0} = sup{i > : \X{t,x)\-i > 0} and Cat, 7 
are constants related to the Wiener process W and respectively to the domain 
O C R\ 

The following notations will be used in the sequel. H = L^{0), p > 1, with the 
norm denoted by | ■ I2 and scalar product (■, ■). H^^{0) is the dual of the Sobolev 
space Hq{0) and is endowed with the scalar product {u,v)-i = (m, (— A)~^w), 
where A is the Laplace operator with domain H^{0) n Hl{0). All processes X = 
X{t) arising here are adapted with respect to the filtration {J-t}- For a Banach 
space E, L^(0, T; E) denotes the space of all adapted processes in LP(0, T; E). We 
shall use standard notation for Sobolev spaces and spaces of integrable functions 
on O. 
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2. The main result 

Definition 2.1. Let x (z H. An H-valued continuous (J- 1)- adapted process X — 
X{t,x) is called a solution to on [0,T] if X e LP{nx{0,T)xO)nL'^{0,T; L^{n, H)), p > 
2, such thatF-a.s. V j £ N, te [0,T], 
(2) 

{X{t,x),ej) = {x,e,)^+J^J^{p\X{s,xmr sign X{s,xm + ^{X{s,xm))Aej{Od^ds 
+ Efcli Aifc Io(^i^' ^)^k, ej)dPk{s), 

For a; e iP(C'), p > 4 and d = 1, 2, 3 there is a unique solution X G L^{0,T; LP{n, H)) 
to ([T|) in the sense of Definition 12. II Moreover, if a; > a.e. in O then X > a.e. 
in X [0,T] X O). 

By the proof of [Sj Theorem 2.2] and [21 Proposition 3.4] we also know that for 
A -> 0, 
(3) 

( Xx^ X strongly both in L^{0,T; L'^{n, L^{0))) and in ^^(O; C([0, T]; iJ)), 
\ weakly in iP(r2 x (0,T) x 0), and weak* in L'^{0,T; LP {ft; LP (O))), 

where X\, A > 0, is the solution to approximating equation 
(4) 

dXxit) - A{^x{Xx{t)) + XXxit) + ^{X^mdt = a{Xx{t))dW{t), 
*a(^a) + XXx + ^{Xx) - on dO, Xx{0, x) ^ x, 

•^x{x) = i (x - (1 + A^o)"^^^)) = *o((l + A*o)"Ha^)), *o(a;) = p\x\" sign x. 

Everywhere in the sequel X = X{t,x) is the solution to ([T]) in the sense of Def- 
inition [23] where X € L'^{0). Below 7 shall denote the minimal constant aris- 
ing in the Sobolev embedding L"+i(C») C H-'^{0) (see ^ below) and C* = 
Sfc°=i Mfc |sfelffi(o) ~ Y^^=i l^k ^k- Theorem 12.21 is the main result of the paper. 

Theorem 2.2. Assume that d — 1, 2, 3 and that 0<a<lif(i=l,2, ^ < a < 
1 if d = 3. Let t := ini{t > : \X{t,x)\^i = 0}. Then we have \X(t,x)\^i = 

0, for t > T, P-a.s.. Furthermore P(r < i) > 1 - (li^^^'y+o ( Jq e-(i-")<^'^ds) ^ . 

In particular, if Ixl^i"' < ^^^^ — , then P(t < 00) > 0, and if C* = 0, then t < 

(l-a)p7i+° ■ 

Remark 1. This result extends to O C M'' with d > 4, if a e !)• However, we 
have to strengthen the assumption on /ife, fc G N, see [Ij Section 4] and in particular 
[6l Remark 2.9(iii)] for a detailed discussion. 



3. Proof of Theorem [272] 

We shall proceed as in the proof of [3, Theorem 4.2]. Consider the solution 
Xx G L^(O,T;L2(f];7ji(0))) to equation (g]). Then by applying the classical Ito 
formula to the real valued semi-martingale \Xx{t)\^i,t e [0, T], and to the function 
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(p^{r) = (r + £2)(i-a)/2^ r e M, we find that 
(5) 

d^e(|^A(t)P_i) + (1 - «)(|XA(i)l-i + e')-(^+"y^X^{t), ^x{X^{t)) + XXxit) + ^xiX^mdt 

_ 1 X-oc ,,2,^ ^^ \Xx(t)ektMXx(t)\l^+e-')-(l~c,f\(X>,(t)e^,X>^(t)).if) 

-■2 2^k = \t^k\^ ^) (|XA(t)|^i+£2)(3 + a)/2 

+ {a{Xy,{t))dW{t),^',{\Xx{t)t,)Xx{t))^, 
{1 ~ a)\Xx{t)el\l^ 



Then letting A ^ 0, by dS]) we get that liminf [ {^x{Xx{t)),Xx{t))dt > p [ \X{t) 

Jo Jo 



Li + °(C')"''' 



-a.s. and hence 
(6) 



+^*/ ([^(g)^^" +^ff(bt)/2 ^^ + 2^ HX{s))dW{s), ^^(|X(s)p_i)X(,s))_i, P-a.s., r<i. 



Next by tlie Sobolev embedding theorem we have 
(7) 

' " i(0)' V u fc ' 'l,L^j, n a > z ana a ^ 



\u\-i < 7\u\lo.+uo), Vug L"+1(0), if d > 2 and a > and V a > 0, if d=l,2. 



Then substituting ([7]) into Q we get 
(8) 

f jJ^^J^^^^W^ds + f {a{X{s))dW{s), ^',{\X{s)W)X{s))^,, P-a.s., r < 



(|X(s)|^i+e2)(i+«)/2 
Now for e ^ we have 



+ (l-a) / (a(X(s))dW^(s),|X(s)|:("+i)x(.s))_i, P-a.s., r <t. 



Hence by Ito's product rule 
(9) 

e-c%i^^)t^X{t)tr + (1 - ")P7'+" r e-^*(^-")n{|x(s)|_,>o}rfs < e-^*(i-")'-|X(r)|Lr 
t 

+ (l-a) / e-c^*(i-")n{|^(,)|_,>o}(a(X(s))dW-(s),|X(s)|:("+^)x(s))_i, P-a.s., r < t. 



From this it immediately follows that e-'^*(i-")*|X(t)|L^", t > 0, is an {Tt)- 

supermartingale, hence = for all t > t. So, ([9]) with r = after taking 

expectation implies that e"^*(i"")''P(T > s)ds < (j^l^^^^, t > 0. This 

implies that P(t X) < ( /q e~'^*(i"")"ds) \ i > 0, and the assertion 
follows. □ 
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